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SUMMARY 

in this paper solutions are presented for the determination 
of transient thermal stress distritutions, in hollew circular 
oylinders of elastic material, in which the temperature is a 
function of the radial coordinate as well as time. Indicated 
solutions are given for an externally insulated cylinder where the 
surface temperature {2 a known function of time, and for the case 
where the temperature-time relations are known for a heating fluid 
within the cylinder. Numerical results have been included fer the 
specific case of a hollow steel cylinder subjected to the flow of 
saturated steam. Yor the particular set of conditions considered 
the tancential and axial stresses were found to be & maximum at 
the heated surface at approximately five seconds after the start 
of heating end attained a magnitude of 52,500 pounds per square 
inch. Stresses of thia magnitude {if often repeated could cause 
fatigue failure and result in rupture of a line carrying any 
appreciable pressure. It is conceivable that for more extreme 
conditions the surface stresses might reach magnitudes high enough 
to cause a permanent set in the material and residual stresses 
would be present when the cylinder was agein cooled to the initial 
conditions. Stresses of these magnitudes warrant consideration 
in the design of steam lines that may be subjected to repeated 


heating cycles. 
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INTHRONUCTI ON 


The thermal strese diatribution in hollow circular cylinders 
censed by & temperature gradient in the radial direction only is 
well known. The case of steady state heat flow in such cylinders 
ig treated in several textbooks. The problem of etrese distri- 
bution for the case in which the temperature is a function of both 
the racial coordinate and time lends itself to solution by the same 
stress ecustions, providing that the temperature-time relations are 
Known, Analytic solutions for trensient temperature distributions 
ere given by many authors for the case where, from an initial uni- 
form temperature, one or both of the surfaces are subjected to 
linear temperature gradients, or are suddenly raised to s given 
value and thus maintained. The process for solution of such prob- 
lems is to substitute into the already evailable stress ecuations, 
expressions for the temperature as a function of time and radius 
for the siven condition of heat flow. It can be seen then that 
for the case of the hollow circular cylinder, as for other confi= 
gureations, the general thermal stress problem is quite straightfor~« 
ward and evolves itself, for the more complex conditions of heat 
flow, into that of the determination of temperature distributions 
as a function of time and distance. 

Consider now the more general problem of the thermal stress 


distribution in a cylinder subjected to the flow of a medium 
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through the cylinder diffcrent in temperature to that of the 
cylinder. In this case the inner and outer surface temperatures 
will vary with time and in general will not be Lineer. The selu- 
tion of such a problem also involves the consideration of a surface 
heat trensfer coefficient which itself may have a magnitude which 
varies with temperature and consequently with time. Analytic selu- 
tions for the transient temperature distribution is a cumbersome 
double series of Bessel functions which does not readily lend itself 
+o numerical results. Hewever, enalytic solutions are available 
fer unidirectional heat flew in the seml~infinite solid in terms 
of exponentiels and the probability function, which are readily 
adopted to numerical results and can be anvlied to the hollow 
circular cylinder within certain limitations, Such golutione are 
valid for infinitely long thin cylinders and give a close approxi-~ 
mation to the temperature distribution for thick cylinders, in which 
the wall thickness is not large compared to the radius. The use~ 
fulness and practicability of such a method to determine transient 
thermal stresses in hollow cylinders when subjected to a heating 
medium is the subject of this paper. 

In the development of the problem, two separate methods of 
attack were considered and thelr results compared. First, a numer» 
ieal method for the transient temperature distribution based on 4 


difference equation rather than the applicable differential equation 
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was used to determine the time-temperature gradient. Seeondly, 

@ purely analytic method wes used to determine the time-tempera-~ 
ture relations and in both instances a numerical integrating process 
was used to evaluate the integrals in the stress equations. It 
might be noted at this point that the stress equatione will give 
exact results, within the accuracy of the numerical integration, 

if the time~texperature relations thus determined are the true 
temperature gradients for the heat flow considered. 

It is the intent of this paper to indicate a practical method 
of solution for such thermal etress problems rather than te obtain 
accurate numerical results for any specific problem. kowever, in 
order to more clearly indicate the feasibility of this approach, 
® simple numerical problems has been carried through, using average 
velues of thermal and olastic properties with slide-rule accuracy. 
The solutions of more complex problems are indicated as the method 
is developed. The problem chosen for its simplicity is that of a 
hypothetical ten-inch I. D. by twelve~ineh 0. D. stemm iine, inie 
tially at 100° suddenly subjected to a saturated steam flow having 
&® temperature of 5007. It is assumed that the steam starts to 
condense immediately and remains at the saturation temperature 
throughout the transient heat transfer precess. The surface heat 
transfer ceefficient was taken ae an average value ever the inner 


surfece of the pipe, considering film tyne condensation, and assuming 
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it not effected by temperature. “he magnitude of the aurface 

heat transfer coefficient was chosen from average experimental 
values for film type condenzation as found in the literature. ‘he 
reason for theee basic assumptions and their validity are further 


digcussed in the body of the text. 
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GENERALIZATION OF PROBLEM AND BASIC ASSUMPTIONS 


The generel differential equation of transient heat flow is 
obtained from the familiar basic Fourier equation for the conduc=- 
tion of heat (4 = ~ AA gt at) and a heat balance. Disrezgarding 
heat sources or sinks and equating the heat gain to the difference 
between the heat anterine and heat leaving an elemental volune, 
the general differential equation expressed in rectangular coerdine 


ates is: 
ey 
BE = ds [2 (eat) + lew gt) + 2 Cena) | 


If k, the thermal conductivity, is considered inverable with 
temperature change and the material is homogeneous and isotropic, 
k will be the same in all directions and constant. 

Then; 


oT = fe eh 2T , 7 
ot CpP oxX* d 4? ee 


& 


In general, engineering materials can be considered isotropic 
end homeseneous and this is an ordinary assumption made in engi- 
neerins calculations. However, the assumption that KE is constant 
may not be valid for all materials within the range of accuracy 


desired in the resulting calculations. Tha order of magnitude cof 
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the variation of k with temmerature may be seen by citing mild 
steel as en evanple, which hee a thermal conductivity of approxi~ 
mately 26 BfU/hr ft °F at 32°F and 23 BTU/hr ft °F et 1100%. 

the desired time-temperature-distance relations for the heat- 
ing or cooling of various shapes are obtained by the integration 
of the general conduction equation, substituting in the necessary 
boundary conditions for the cage in question. One of the simplest 
cases is the unidirectional heet conduction in a plete of semi- 
infinite thickness, with the free surface at x xe. If the 
term k/eye is defined as the thermal diffusivity and indicated by 
the symbol X the general conduction equation for unidirectional 


heat flew reduces to: 


oT — Awe 
ot = (4) 


fo solve this ecuation, it is necessary to determine the 
formulae which will express the initial end boundary conditions 
which the temperature satisfies. If a semi-infinite solid with 
ite plene face at x = 9, initially at the constant tempsrature of 
sere end suddenly subjected to a heating medium, which has allied 
with it ea surface hezt transfer coefficient E, the initial and 
boundary conditions ere: 


Initial condition "T =e WHen t=0 


Boundary ecendition out 2. (Ts - Ta) = © AT x=0O 
ax Je 
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or defining H/k as h the boundary condition is: 
ane vt + (Ts - To. ) =O Aer xX=O0 
x 


The solution of the differentia! equation for unidirectional 
heat conduction with these initial and boundary conditions will give 
the temperature es a function of time and distance for the semi- 
infinite solid. Such a solution will alse be valid for a plate 
of finite thickness having a large ratio of surface to thickness, 
if it is heated, as before, at the face x = @ and the other face 
is perfectly insulated. To use such a solution for a hollow cir= 
enular cylinder, externally insulated, and subjected toe the flow 
of a heating medium it must be as@umed that the curved section 
of the cylinder can under certain circumstances be considered as 
a flat insvlated plate. Thies asaumtion will necessarily introduce 
an error in the transient temperature distribution thus determined 
for the cylinder and will be a function of the ratio of thickness 
to redius. This fact muct be considered in the caleulation of 
stresa diatributions for stresses thus determined are only as reli~ 
adle as the temperature distributions used in their caleulation. 

To verify the validity of this assumption the applicable 
difference equations may be set up for the plane plate and the 
hollow circular cylinder and the results compared. As explained 


under the heading "Numerical Solution of Temperature Distribution®, 
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the numerical method converges rapidly to the analytic solution. 
this being the case, a comparison of numerical methoda should alse 
justify the use of the analytic aclution for the plane seni-infinite 
solid for stress calculations for the hollow eireular cylinders 
within the same range of accuracy. | 

Taking a section of the plane solid of unit width, Ax thick, 
and dividing it wp es shown in the figure below, & heat balance 


gives: 





FIGURE A 


k(T.-Ta)at + k (T3-Ta) at = CoO Ar? (Tz - Ta) 
WHERE: To) = T2 Ar (t+at) 
Oe (T-h) + (T3-T2) = ™ (T-T2) 


WrHEeee: M = CpPax* 
Hw AT 


ny 
j 


THEN: = T2+ h -2T + Ts () 


wa | 
T + G-2%H+13 Wee M=4 
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A similar heat balence for the cireular section where the 


ares betveen sections progressively increases gives: 





FIGUEE & 


Ae (T.-R)at + & (Lo + avae)(T3-h) at = CpPaxt(t + ax00)(Tor -T2) 
WHERE! Toy = la Arter (t+at) In CyLinDeR 


OR (T-Ta) + (10 +axae)(Ta-f) = M (1 + 4x40) (Ter - Tz) 
Z 


SOWING For ee = LETTING Me 4. 
ie em + \n- - T2 eee + Ts (1+ avao) (2 
(\ ea 


It can bs seen that equations 1 and 2 are identical in forn, 
the only difference being the addition of the term 4x40 for the 
cireular cylinder. For the ordinary problem, in order to satisfy 
the basic assumptions of such & numerical solution, Ax mist be 
chosen as & low order term. This means that the radius of the 


cylinder will be large as compared to Ax and thus Ae must also 
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be a low order term. The magnitude af these values may readily 
be seen by taking as an example a hollow cylinder of five inch 
inside radius and a ax = .01 ft (Note that this is the largest 


value of'ax chosen for numerical calculation) then: 
A® = TANAO = QA? = .024% RAD. 


and Ax4e 2 .00024 (Yhich is essentially a 


order term). 





Then for the use of a modulus M 2 4: 


( 
ac sc + E - (2.00024)T + (1.00024) Ta | 8B) 
4.00048 


A comparisen of equations 3 and 1 shows that this point, 


namely at a radius of five inches, TS and tT) 


identical. However, in the numerical solution for temperature 


. are practically 


distribution, the temperature at any point is used to determine 
the temperature at the next adjasent point, ete., through the 
solid. This means that the error in using equation 1 instead of 
equation 2 will be an accumulative error and places a restriction 
upen the thickness of the cylinder for a given error. This fact 
perhaps can be seen more clearly by noting a Figure A thet; 
AXn =AX+ HAxXo6 
Wheres n is the number of steps te be taken 


for the given thickness chosen. 
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It should also be noted that the error is not a function 
of Ax alone but a function of the product of Ax and Ao, 
which essentially meane that for a civen radius cylinder and wall 
thickness the error involved will be the same regardlecse of the 
choice of ax. With these facts in mind, it 18 apparent that the 
limiting factor is the ratio of wall thickness to radius regardless 
ef their absolute magnitudes. 

Another pertinent observation is that at the first instant 
of heating the two equations almost coincide and the error in 
temperature distributions increases with each step. As the temera~ 
ture distributions —— determined, only enter the streas squations 
in the terns [trae and Trolr the true criterion for 
the error tevetved in atreas tintiite would be to plot tempera- 
ture times radius vs. radius eurves for both methods and compare 
their areas. 

It 4a the author's opinion that the error involved due te this 
vasic assumption will be of the same order of magnitude or leas 
than the errors inherently involved in thermal and elastic proper= 
ties, where ratios of wall thickness to radiue found in ordinary 


atean lines are to be ueed. 
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WYMERICAL SOLUTION FOR TEMPERATUME DISTRIBUTION 


Unidirectional heat flow will eccur in a flat slab at a 
distance from the edges and alse in an insulated rod of uniforn 
gection. As far as temperature distribution is concerned, the size 
and shape of the rod ie immaterial and ia only of consequence in 
calculating total heat flow. Consider such & red of unit width, 

a thickness Ax and divided wp inte sections Ax long as indicated 


in the figure below: 





If now a finite surface heat transfer coefficient H is con= 
sidered between the surface (0) and the ambient fluid temperature 


(a2) then for the point (0) a heat balance gives: 


ri ax (Ta ~ To) Aca = k (T —To at + CoP ax? (ie 1a) 


Wheres To! =< To Ar (t+at) 
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The finite difference equation thus obtained is based 
upon the following assumptions: 

a. That At be taken amall enough so that there is a neg= 
ligible error in heat flow calculated by the tern 
(7, -To)/axe 

bd. That At be taken small enough so that there is a neg 
ligible error in neglecting any but adjacent areas 
(ax by ax). 

ec. Thatax be taken small enough so that the tempersture taken 
at the center of the element can be used to commute the 


heat capacity. 


Define: ME nGele ans 5 Bie 
RAT K At 
N= Har 


Using these relations: 
ive ts) + (het) = m4 (ge ~ To) 


Solving for T} 
we Sooner ¢ T; 
o = 2N +(\- 2N4 + 
ma ante) : 2 
This equation may be written in terms of "influence factors" 


as follows: , : 
US Paola + lee To + eae 


Where: 
ie) a ) 0° é a } } O\ o— 
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The subscripts indicate in order, the influence that the 
present temperature has on the point whose temperature et (t+ At) 
is to be found for the period of time At. In other words, if the 
new temperature at point (0) is te be found after the time At, 
Fag indicates the influence of point (a) on the temperature at 
peint (0) during the time interval At. 

If f3 is to depend upon f, in a positive sense, then: 


Ess = |\~ ONG 2o 
~ 


Or M-2N+2 26 
mis 242 

Since all of the quantities entering into KN are inherently 
positive, it is apparent thet M met be larger than two , for the 
case where a surface transfer coefficient is involved. Depending 
woon the magnitude of the surface transfer coefficient and the ther- 
mal constants of the material, this may place a limitation upon a 
convenient averaging process. The value of M=4} was found suitable 
for numerical temperature calculations for steel, 

The general mothod ia to use the weighting process, as described 
sbove, for surface temperature caiculations and the simpler avrera# 
ging process for interior temperature determinations. It igs con- 


venient to choose a scale such that fractions can be neglected. 
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If the fraction 1/2 was consistently added or dropped, an error 
would tend te aseumilate. Therefore, a rule of dropping the 1/2 
efter even dicits and adding 1/2 after odd digits can be adopted, 
ass 1/2<0, 3/2-2, §/2-2, 7/2-%, ete. 

The choice of the modulue M is diseussed in reference (5). 

The larger the modulus M the sooner the numerical results will 
converge to the analytic solution (4.6. in leas steps) and the 
greater the mumber of calculations necessary to cover the same time 
period. After due consideration of the work involved, & modulue 
M=4h was chosen to obtain rapid convergence to the analytic solution 
during the early time intervals. 

Now consider the heating of a one~inch thick plate insulated 
on one surface and subjected to a heating medium on the other sure 
face. For simplified calculation, censider EK and H ag constant. 
By choosing several values of ax (and the compatible at) for a 
modulus M=4, surface temperature ve. time curves can be plotted to 
determine the validity of the choice of the magnitude of az, As 
the accuracy of the solution is based on the assumption that ax 
be emall, such a plot should indicate any variation in the results , 
due to an increase in the magnitude of ax. 

For ordinary Steel: 
2.50 f¢"/ne 


k 225 Bru/hr, ft, F 
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Assume: H= 1000 Bru/ hr, ft”, So (A fair aporoxima- 
tion for film-type condensing steam ~~ set Appendix 1) 


Semple determination of ¥, Fao? Feo Fol and At. 








a? 2 — vowrT. AX 712 x.01 = 112 INCHES 
Then; N= Hax = 1t000%x%.01 = .40 
Ae 25 
M= ax® Om et et WHere M=4 
Kat 4K 
= (.01)* = 6.5 xIo7* HRs. 
4x.5 
At = 6.5% 10°" s Sooo = 18 SEC. 
~ 4 
feo = (1 - Ziti) = (> aibieaio) -o%8 
m 4 
es = Z = a - 0.85 


i 
Similarly for the other values of AX: 
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The eqnations for the surface temperature become; 


— soem + SSem% + 087, 
™~ = Siem + ofS Tr + o8T 
a Tas olom™ + eet + ost 
™ Te = 005TH 4+ 045 To + 0.57, 


On set of calculations will be given at this point to indie- 


cate the form and precedure. The results only for other numerical 


calculations are included in Appendix II and pletted on the enclosed 


curve sheets. 
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ANALYTIC SOLUTION POR THHPERATURE DISTRIBUTION 


fo determine the time-temperature-—distancs relations in the 
gemi~-~infinite solid for the case where the initial temperature is 
zero and the surface temperatures is a function of time (+), the 


following equation is valid: 


of 
p2 
T = 2 b(t - 2 je da WO 
tr 4K pt 
X/2VKt 
Also for the caso of the semi-infinite solid, initially at 


sero temperature, which is heated at the surface x=0 by radiation 


from a medium at constant temperature T,, the applicable equation 


ie: 
hx + h* Kt 
" jeri Le eric / 2% . VRE (2) 
Ta 2VKt 2VKt 


Where erfe is defined ae one minus the probability function 


or error integral. 


evfex = ( i— 2. 


ale 
Je) 
te 
0 
3 
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It can readily be shown by direct substitution into the 
diffarential equation of unidirectional heat flow that theese 
equations do setisfy the differential equation and initial and 
boundary conditions. 

Equation (1) was civen to emphasize the rr that whenever 
the gurface temmerature of the gubetinan pe wien as e function 
of time $(t), and this function is continuous and can be integrated 
ever the time interval (in equation (1)) to be considered, the 
resalting time-temperature-distance relations for the solid are 
then determined, This opens the avenue for the solution of a number 
of interesting problemas. As noted before, a numerical solution 
for surface temperature is readily adopted to thermal provertices 
of the material which vary with temperature, as well as a surface 
heat transfer thet is a function of temperature. This method can 
also be ervtended to the case where the temperature of the heating 
nedium is some known function of time. however, it should be noted 
that such a eolution for surface temperature, under the above stated 
conditions may become very tedious, as the influence coefficients 
may have to be calculated for each atep of the procesa as well as 
choosing a At that ia small. Althongh the numerical solution 
iteelf, if carried through the solid to the desired depth, will 
give the time~temperature-distance relations directly, the desired 


surface temperature $({t) may be determined over a shorter interval 


with the resulting saving of much numerical work. 
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Figures (C) and (D) show the determination of the surfaos 
temperature relations O(t) for a short interval of time, Figure 
(D) shove that the surface temperature ve. time curves very closely 
approximate a straight line on a log-log plat, and can readily be 
pat inte an anelytie expression. 

Equation (2) has been used to determine the time~-tomperature= 
digtanee relations in a perfectly-insulated, one-inch-thick plane 
solid of lerge surface area. The resulte of thie can be used as 
a close approximation of the true transient tempereture gradient 
in an insulated LO-ineh inside dismeter hollow cylinder of one~inck 
wall thickness. 

For the ultimate caleuletion of the stress distribution, it 
te convenient to plot tesmerature<distance relations for various 
intervals of time, in order te perform a graphic or numerical inte- 
eration of the integrals in the stress equation. Figure (8) shows 
a plot of these relations and shows how closely the numerical method 

arallele the analytic solution. Figure (€) gives the equivalent 
temperature times distance vs. distance relations, und indicates a 
graphic integrating vrocess, It might be preferable and possibly 
& little more accurate te use a numerical integrating scheme esush 
ag the travezoldal or Simpson's rulis. 

Yumerical resulte ore not as tedious as a first inspection of 


equation (2) might indicate. Tabulated values of the error integral 
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are available and many of the parameters are eonstant for one 
Snetant of time. A simple tabulation of the quantities given 
below will materially facilitate calculation, all other para~ 
meters being constant for a given inatent of time. 


Gi x wr 








(2) x 
Vt 
(3) ey fc Wee 
(4) evfc (Re + h Vee) 
(5) hx 


(6) hx + h*Kt 
(7) e+ nikt 


(8) : 
The reanlts of calculations made for the semi-infinite solid 


for a penetration of one inch are given in appendix te correlates 


the data used to plot curves. 
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DETREMIMAPION OF STHESS DISTHIBUTION 


If the transient temperature distributions determined by either 
numerical or analytic methods for the semi-infinite solid under the 
conditions specified are considered good approximations of trues 
temperature distritutions in the hollow cylinder, calculation of 
stress distributions are atraight forward. At any instant of time 
the stress becomes a function of the thermal and elastic proverties 
of the news for the ceometric configuration considered, and can 
be readily calculated. The following thermal stress equations for 


the finite hollow cylinder are given by Timoshenko: 


Keni [rw - [rae | 
brat 
Je = ; : 4. ae _ ts — Tr 
as +. rte Trdr - om 

(=¥) ©) meas 


These equations are based on ue and = properties 
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that do not vary with temperature. Within this linitetion they 
are exact for an infinitely long cylinder of any wall thickness 
and any given temperature distribution so long as it is in the redial 


direction only. The principal planes of atreas are in the radial, 
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axial, end tangential directions with normal stress acting on 
each and stress acting along a plane perpendicular to this axis 
are independent of the position of this plans. 

In order that the solution be exact when a cylinder of finite 
length be chosen, boundary conditions require that external forces 
equal in value and distribution to the stresses in the axial direction 
act on the ends. In the ueval case such forces will net be present, 
eo that the stresses near the ends will be very different than those 
calenlated by this methed. This effect diminishes eas the distance 
from the ends is increased and practically disappears at a point 
sufficiently far away. Writere differ as te the distance from the 
ends where this point is lecated and it actually depends upon the 
ratio of the thickness to cylinder diameter. For the purpose of 
our problem, this influence will be considered negligible. 

Ig is interesting to note that in the case where the surface 
of a solid, initially at uniform temperature, is suddenly raised to 
a given temperature and thus saintained, the maximum etress occurs 
at the surface at that instant. In fact, the meximam stress (i.e. 
tangential and axiel stresses which are equal at the surface) will 
always occur at the heated surface for a solid initfally at uniform 
temperature. The stresses thus produced are equivalent to the stress 
produced in a thin plate for the same temperature rise in which a 


chance of area is prevented. This stress is only momentary and 


Zit 


rapidly drops as the increase of temperature spreads inward and 
allows an expansion of the outer layer. Most published solutions 
have this common boundary condition, namely, an instantaneous 
surface heat source. 

For the case of radiation into a solid from a fluid, where 
the surface temperature is a function of time ard equal to the 
initial solid temperature at the first instant of heating, the 
conititions are vastly different. The maximum surface stresses 
will not occur at the first instant of heating but at some later 
time, and depending upon the conditions of heat flow may even occur 
at the steady atate condition. For the case of the sudden rise of 
the surface temperature of a solid, initially at conetant temera- 
ture and ultimately coming to a steady state condition, it is essen-~ 
tially going from an unstressed condition to a streseed condition. 
In the problem chosen, however, the cylinder at an initial uniform 
temperature with the outer surface perfectly insulated will, upon 
heating of the inner surface, be taken from one unstressed condi~ 
tion to another as the cylinder finally reaches the temp*rature 
of the heating medium. It is reasonable to essume then that maxi- 
mam atrese conditions will occur at some intermediate time, 

Caleulated values of inner surfece strees are tabulated in 


Appendix If for a one-inch thick, ten-ineh I. D. circular steel 
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cylinder along with the strese distribution across ite section at 
the time of maximum stress. The cylinder being insulated on its 
outer surface, initially at 100°F and auddenly subjected to a satur= 
ated steam flow which is maintained at 500°F. The results of these 
esleulations are graphically represented by figures (£), (~) and (G) 


and will be further diseuseed in the results and conclusions. 
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APPEHEDIX I 


Acenrate thermal streas determinations in eteam lines are 
usually linked with a knowledge of the surfece heat transfer co~ 
efficient. The matter of the vhysical phenomens involved as well 
as the magnitude of sucl/ coefficients is the subject of a good 
deal of conjecture. This coefficient is defined by the heat flow 
equation: 

q = HA(T,-Ta) 
Wheres q - Heat Flow 
HK - Surface transfer coefficient 
A ~ Transfer area 
Tg~ Fluid temperature 
Tg- Solid surface temperatura 

Values of H determined by experimenta] means are in general 
more reliable than those determined analytically, and particularly 
for condensing stean. 

there are two types of condenration, namely, dropwise and 
film condensation. Either or a combination of the two types may 
exist, depending largely woon the cleanliness of the fluid, clean- 
liness of the condensing surface and the surface conditions of the 
material concerned in the condensation process, In any case, Leat 
transfer coefPleiente for the dropwise presess are in general many 
times higher than for film condsnaatian, but are more difficult 


te determine, 
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Authora differ widely as to the magnitude of the difference between 
the two types, giving coefficients from 4 to 20 or more times larger 
for the more efficient dropwise precess, It bas been determined 
that in ceoneral clean «team with er without non~condensidle gases 
always condenses in a film en clean surfaces, whather rough oF 
emcoth, 

Experimentally determined valueea of heat tranafer coefficiente 
are given in terms of an averase value over a prescribded ares, 
discounting the fact that the value of the coefficient may vary 
considerably over the area depending upen the film thickness. Nusselt 
in an enalytic investigation of heat transfer coefficients on eyl~ 
indrieal surfaces found the film thickness and consequently the 
trensfer coefficient to vary arcund 4 horizontal tube, His analysis 
shows the mininwm film thiekness and consequent maximum transfer 
coefficient to exist a few degrees on sither side of the top of a 
horigontal tube. Kis investigation ia equally valid for condenga- 
tion inside or ontside of a tube. By integration around the tuds 
an average ceefficient was determined, Experimental determinations 
heve shown the average coefficient thus determined to be somewhat 
in error and correction factor curves are available for use with 
Nusseli's equation to give representative values of the surface noat 


transfer coefficient for condensing atean on tubes. 
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If now dropwise condensation were to be considered, the 
stress problem imaedinately becomes very involved for tha lecal 
tranefer coefficient would vary avvreciably from droplet to un- 
wetted surface and the heat flow would se longer be uniédtrectional. 
For auch a problem an entirely different approach would have to be 
considered. ‘The value of H-1000 BYru/hr ft” OF used in the numerical 
ealculations is an arbitrary value chosen from the references for 


film-type condensation of steam on & horizontal tube. 
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RESULTS AWD CONCT IIS TIONS 

It is apparent fron the magnitude ef the inner surface strese 
determined, for the numerical problem considere€, that thermal 
siresses uncer severe Sonditions of repeated heating eveles could 
contribute to ultimate rupture of a line under pressure. The tan» 
sential atresa on the inner surface due te preagure in the Line 
will be tenelle and to some degree counteract thermal stresses. 
However, for the 10-inch I, D. by l2-insh 0. DL. line considered, 
the stress due to a pressure of 515 #/in® in the line would be ap~ 
proxinately dr = 500 ¢/ in” in compression and do = 1,6iW #/ in* in 
tension, These atresees are extremely amall in magnitude as com 
pared to the thermal stresses determined. 

it should be remembered thet this peper was motivated vy @ desire 
ts formulate a practical approach to the actual thermal stroes 
problems in hollew cirenlar cylinders which are extremely difficult 
if not impossible by purely analytic appreach. Recesserily, many 
esgumptions were made, and to some degrees they have been justified 
in the text of the pamer, It is obvious then that the rewults of 
the numerical calenlationag should be substantiated by an akper Laveane 
tal investigation before they are eecepted as accurate enough for 
practical uege. 

Aas little data 1s available on the magnitude er variation ef 


surface transfer coefficients, fer condensing stom in hollew 
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cirmular cylinders under high temperatures and pressures, the 
value chosen covld be considerably in error, as it is a value 
extrapolated from experimental ree:lts for condenser tubes under 
Lower temperatures and pressures. It is conceivable that a 
consideration of the ehance of the magnitude of the surface trene- 
fer coefficient with temmerature difference across the film as 
well as the change in thermal properties with solid temperature 
would materially effect the thermal streeses determined. The 
simplifying aesumotion that the line was perfectly insulated would 
aleo constitute some error in stress caleulations beyond the tine 
at which the outer surface temperature begins to rise. Although 
uaier steady state conditions ingulations are available that are 
extremely effective, under transient conditions, temperature dise- 
tributions will necessarily be effected hy the heat capacity ef 
the inguletion, This condition wes net considered in the nunerie 
eal problem which injects sone error in the reauiting etreas deter~ 
minations for times beyond approximately five seconds. Perhaps, 
with some of the qualifying assumptions in wind, thia paper will 
atimilete further investigations along the zamm linea toward a 


useable practical galution for thermal stresses in steam Lines. 
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RRO OMMENDATIONS FOR FURPHER SIUDY 


Further study of the practical problem of thermal stresses 
in hollow eirculer cylinders during actual tranalent boundary ecne 
ditions falls into the several broed categories listed below: 

(a) The experimental determination of ths true boundary 

conditions thet may exist in the actual problem during the 

transient state, such as these conditions that may exist 
when a eateam line 4s suddezly subjected to a flew of high 
temperature stean. 

(>) An experimentel verification of the numerical example 

given in the text as to the time at which maxinnm stress 

ocours, its magnitude, as well an the streas distribution 
across the section at time of maximin streas, 

(ce) An experimental study of the Magnitude of the surface 

transfer coefficient under high temperatures and preesures 

in hollow circular cylindors. How the value of such a co~ 
efficient varies with presaure, temerature difference 
across the condensate film, and anezlar vositian in a horl- 
zontal eylinder. 

(4) A further fustification of the use of temperature dis- 

tributions as determined for unidirectional heat Plow in 


the stress calculations for hollow cireuler cylinders, 
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There are two immedints averuca of approach to the last 
category, firat, by direat comnariaorn of the teaperatura distri- 
butions as determined by purely mamerical methods fer the plane 
plete end hellew circular cylinéer, and, eecondiy, such & tempers 
ture comparison could be made by the use of a conformal mapping 
process from a plane finite section to a circuler section. A plot 
could then be made of the error in stress calevlations as a funce 
tion of the ratio of wall thickness to radius, which would verify 


the usefulness of the approximate method for e problem of given 


geometry. 


30 







Mant wil of er _ — ah - 
ahthth omberepet wit © 0 : te 
| =e ret 
~m_s * ae = my 
oo y= y= oa AL 
guy hve haves salerais pou 3 die! 
Mesto ke winnie stent wt in 


Tilwe bbws os ausmatads 
ete & abideng 2 20? 0? Bedtime eat iver 


- ©& 





* 
~ 


























Ya ae 





—— Le LST | 
- ——&—— - « 





sée -_ —_—<—- ——-oa 
7 - ms — = 


a ame ol 


: ia* ("am ~ 
— ———-_ =P 
' »« =e ee 











+ —— >: O00O0O8Oereww a 
| Sueeace TERE RATURE 

———ae | her o 

TIME 
3h eet, ere 


~—T — 


—_— ae . > * a . 


_ Ph De} — 












= - 
- &. 
. 





€ fer 


= — if 1 | : 
_ Se j _ | 


Mere TEM ATV? 





- & 
© (2) 9 
a Q)-~ 06 “” ove" 
~ OF 4 Z§- 
v A> ig Ne a 25 
ee —_=—— 


ro. 







"aos WI 


























20L10 TEMPERATURE 






JIBO 






NurqtéreicAd Ano ANALY TIC 








IVETHODS 


Jet 


—— 


ANALYTIC DoeuTioNn (Socio Cveves) 


JZ =erf =e sees ah fs md 


7a, 2’xt 


NUMERICAL PIETHOD (PAsHED Coemea) | 


Suerace Ten Je = 02k +03h +057 sd 


Sewo Tene EZ =7+24a+kh 
- 


a a 


—EE —. 


RAOIVS IN INCHES 


aL <1 9/4 a 


Terriat Tera -100°° Hearinc (Meépie -Scecr =§9H=1000 


- 
— ‘ —_ 2 —_ 













SURFACE JEMPERATURE (no A= 


Ny 


{20 


SUELACE TEMPERATURE 





Vs > 
TIME 
= I —- - ——E 
| 
mn ante 2 
; 
aaa s : 
— — i 
= —_—_—+ . = : 





| 
; x 
a. = @ 


’ 
| | 
—_ ~eecy Al se | 
0 20 Os OC SO 70. OD, BD tO 


TIME Iv SECONDS 


FICURE C Devt -InFIniTe DeLiD 
fmrinag TevtR-100%  HeATING (PIE WSrI +500 


— 


ei 





| =) | 


| Suerace Taewwenrune Vs Tae i 
GL, ha 


| 
| 





ar o2 of 26 ar a8 OF Ne 
—- is snc aves 


FIGURE D Sert-inFunirTe Sowioe 
LmiTiAd TErtR -100°F Hearing Meowrt ~ 500 °F 








a 


TEMPERATURE Ties RADIUS 


(200 - aay Vs = 
\ RADIUS 

/200- \ : ) , . ; 

Ms \ 

Sas 

yw 


8 






= : Bs 
S02 SA 5.2 53 54 SS 54 5.7 5.8 45.9 &0 
RADIVS /M INCHES - 


Frovurée F jo" id Br 12"00 C Yer OER = IniTrAis.y 
AT 100°F SVOPENLY BYBJECTED To THE Fiow 
OF S00°F SATURATEL STEAL. 










FiGur& F 10" 1D 6r (2°OD CYLINDER 
ArT 100°F SevoPen.y SuRBSs C ED 
or So00°F DATYRATEP OTEAC. 


Se, 000 - + 
| | [a Fibs 
ee | 


STRESS DISTRIBUTION 
ATTINEORF 
Maximum Ireess 






ConPeessive& Sreess - */i? 


TeNsate Sreess - thn? 











a a 


WURLO GIIVIOLES 700s YO Mors FH of a sek So 
AINFIICS F.00F LY APTWHIINT — FIANITAD J0y2) 4907,,01 2 43909). 
< | = 2 
| SONODIS V/h Phidf se any 
eof 68S OR OK, BF, OF CR _ CR eee 


r as me | = 





— =) 








-_— »@ — 
( ' s : 
| ’ 
— —— @ *. —— st 
“ ak 
’ ; 


tk " DNL SA csctaic sorvve27e 


| | | ye emt te 





APPENDIX If 


Tabulated values used to plet curves 
Figure (A) — (Numerical method ~ Surface temp. ) 


Tez Ta + (NETS + 4% (Nore: =A) 
> pa 


Te=!00, Ta = 500, W=tooo, K=.50 


Ax = .12 inches ax = .09 inches Ox = .06 inches Ax = .03 inches 
at * .180 sar At = .101 see At = .045 ses At = .01125 sec 
Tin’ ¢ in Seo in®F tin See Tin? + in Ses Tin°F % in See 


100 “9 106 0 100 0 100 0 
120 18 160 101 140 0 O45 129 0.01325 
204 036 181 202 156 0096 129 0.02250 
221 o kS 196 » 303 167 0135 137 0.03375 
234 272 208 o lvtt 178 o L86 142 0.045 
245 09 217 0 5OS 185 oLe5 | 146 0.05625 
226 » 606 192 0270 159 0.06750 
233 2707 198 2315 154 0.07875 
2h0 - 808 203 «360 
216 909 208 «405 
212 450 
216 2 4O5 
Z20 540 
2ah o 584 
227 » 690 
250 >675 
£33 » 720 
2°36 765 
239 e810 
ahe O55 
ahip » 700 
27 0 G45 


Figure (A) - (Analytic method ~ Surface temp. ) 
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To = 100%, %, = 500%, h = 40, K = .50 


\ 


f in tin See 
143.0 0.05 
159.6 0.1 
177.6 0.2 
192.4 0.3 
204.3 0.4 
215.2 0.5 
221.7 0.6 
229.7 0.7 
236.4 0.8 
242.0 0.9 
247.0 1.0 


Figure (8) ~ (Humerical method - Solid temp. ) 


jetties tap. = To = O2Ta + 0.387, + OST 
Selid tem. — Ta = T s 272 +713 
A 


f, = 100%, tT, = 500%, H = 1000, K = .50 
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50.9 See t2l.8 Sec t82.7 sec tH2.6 See tae Seg 


Radius 
met , - 8 f — 8 tr-% # £wt-% tT - YF 
0 Zh5 281 303 318 331 
012 161 203 231 251 26? 
» 2h 117 151 177 198 215 
236 103 121 My) 159 174% 
2448 100 107 119 133 145 
60 102 108 117 127 
072 100 103 107 Lis 
o Sh 101 103 109 
96 100 101 103 
1.08 100 101 
1.20 100 
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Figure (@) - Analytic Method ~ Solid Temp. 


t = 0.1 See 
Radius 

T- % in * 
158 5.0 

136.9 5.02 
124.6 5.04 
1144.1 5.06 
107.6 5.08 
103.8 5,10 
101.8 5.12 
100.8 5.14 
100.26 5.16 


7, = 100%, f, = 500°F, H = 1000, k = 25, h™ MO, K = 50 


t = 0.2 Ses 
Radius 
g.- °F in ® 
174.8 §.00 
146.2 §.03 
130.7 5.06 
116.4 5.09 
108.7 5.12 
104.0 5.15 
101.7 5.18 
100.7 5.21 
100,39 5.24 


t = 0.5 See 
Rediue 

t-~ FF in * 
214.4% 5.00 
188.8 5.03 
166.2 5.06 
148.0 §,09 
132.5 5612 
122.9 5245 
115-2 5.18 
109.6 5.21 
105.8 5.2% 
103.5 5027 
101.9 5230 
101.0 5.33 
100.5 5.36 


t = 1.0 See 
™ Radius 
fT~o F in* 
247.5 5.00 
208.7 5.05 
176.8 5.10 
152.0 5.15 
134.0 5.20 
126.5 5.25 
112.1 5.30 
106.8 5.35 
103.6 5.460 
101.7 50'S 
101.1 5.56 
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4 = 3.0 See t @ 5.0 See t= 7.0 Sea t = 20.0 See 
Radius sedius Readius Radine 
T-~ %F in ® tm yy tn ® p . %F in ® g++ % in ® 
307.5 5200 335 5.00 354 5,00 370.5 5.00 
248.5 5.10 462.5 5510 307 5010 329.5 5.10 
199.9 5.20 237.3 520 266.3 520 293.7 520 
163.5 5.36 195-5 5.239 227 0F 530 2552.90 5.32 
137.7 540 169.8 5.40 197.0 5ei 225.6 45.40 
120.8 5259 *h6,7 5.50 170.0 5.50 198.8 $50 
110.9 5.60 130.4 5.60 151.2 5.8 17509 5.8 
105.3 5.70 119.3 5.70 13%.7 579 157.6 5.70 
101.7 5.80 110.9 5.80 124.9 5.80 WB 5.80 
100.7 5.90 105.4 5.90 125.1 5.90 192.0 45.90 
100.2 6.00 103.3 6.00 107.8 6.00 124.5 690 





Figares (C) and (D) ~ (Numerical Method - Surface Tem.) 


To' = NT + (2 -NE)TO + tT. (Nore m=4) 
AX = .09, at = 1ol , Ts = loo, Ta = 500, Ge".50 

Time in 

Sec H = 900 B = 1000 _= 1100 H = 1200 

0 100 100 100 100 

101 154 160 166 172 

202 174 181 1&8 195 

» 303 188 196 204 211 

hoe 199 208 216 224 

0 505 208 217 226 235 

- 606 216 226 235 2h 

0707 223 233 242 252 

2808 230 pt ale a9 259 

909 235 246 255 265 





Fieure (E) ~ (Analytic solution « Temperature times radius 


ve. radius) 
fadius t = 1.0 See ¢ = 3.0 Sec t = 5.0 See t+ = 7.0 See t = 10.9 Sac 
in 

oe = a bee bee... BAS. Lhe 

733.0 1038 1175 1276 1352 

3492.0 758 931 1055 1279 

177.9 520 Ty 855 1008 

64.0 337 527 676 B21 

19.5 293 377 524 675 

6.5 114 257 385 eral 

61.0 170 287 125 

30.2 120 198 328 

9.8 58.0 144 260 

4d 31.9 B9.0 188 
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Figure (G) ~ (Analytic Solution - Surface stress vs. tine) 
Surface 
stroas t=1.0 Sen $23.0See $25.0 See t57.0 See  t=10.0 See 
in 
Los/in® ~-38,950 wly9,700 50,250 aly , 800 v7 ,800 


Fisure (F) ~ (Analytic Solution ~- Stress distribution at time (5.0 sec) 


of maximum stress) 


in Je or Oz. 
Inchea_ Laa/in® Loslin® Lbs/in® 
200 $0 250 0.0 50,250 
5.10 33,700 B22 34 9522 
5.20 19,639 1,192 20 ,822 
5430 6,510 1,530 8,340 
540 B26/ 1,573 949 
5-50 7,650 151 6,199 
5.40 1Z,330 1,243 11,887 
5270 15.550 965 Lr, 585 
5080 17,900 659 17,242 
5.90 18,900 333 18,567 
6.90 19,216 9.9 19,219 
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